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Abstract 

In this paper, we extend the calculation of tensor vacuum susceptibility in the rainbow-ladder 
approximation of the Dyson-Schwinger (DS) approach in [Y.M.Shi, K.P.Wu, W.M.Sun, H.S.Zong, 
J.L.Ping, Phys. Lett. B 639, 248 (2006)] to that of employing the Ball-Chiu (EC) vertex. The 
dressing effect of the quark-gluon vertex on the tensor vacuum susceptibility is investigated. Our 
results show that compared with its rainbow-ladder approximation value, the tensor vacuum sus- 
ceptibility obtained in the BC vertex approximation is reduced by about 10%. This shows that 
the dressing effect of the quark-gluon vertex is not large in the calculation of the tensor vacuum 
susceptibility in the DS approach. 

PACS numbers: 12.38.Aw, 11.30.Rd, 12.38.Lg, 24.85. +p 
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The QCD vacuum susceptibilities play an important role in characterizing the non- 
perturbative aspects of QCD and in the determination of hadron properties M- In par- 
ticular, tensor vacuum susceptibility is relevant for the determination of the tensor charge 
of the nucleon in the QCD sum rule approach [4]. The previous calculations of the tensor 
vacuum susceptibility have shown that the theoretical treatment of this quantity is subtle 
and different treatments can lead to different results . In order to get a reliable theoret- 
ical prediction of the tensor charge, one needs to determine the tensor vacuum susceptibility 
as precisely as possible. Recently, a particular implementation of the vacuum polarization 



definition of the vector vacuum susceptibility has been proposed in Ref. [10|, in which the 
vector vacuum susceptibility is totally determined by the dressed quark propagator and the 
dressed vector vertex. Soon this definition of vector vacuum susceptibility has been general- 
ized to calculate the tensor vacuum susceptibility by some of the same authors in Ref. 
Just as was shown in Ref. ^|, in order to calculate the tensor vacuum susceptibility, one 
needs to know the dressed quark propagator and the tensor vertex in advance. At present 
it is impossible to solve for the dressed quark propagator and the tensor vertex from first 
principles of QCD. So one has to resort to various nonperturbative QCD models. In the 
past few years, considerable progress has been made in the framework of the rainbow-ladder 
approximation of the Dyson-Schwinger (DS) approach 12|-|l7l|. Due to the great success 
of the rainbow-ladder approximation of the DS approach in hadron physics, the authors in 
Ref. jlll adopt this approximation to solve for the dressed quark propagator and the tensor 
vertex and from these obtain the numerical value of the tensor vacuum susceptibility. How- 
ever, it is well known that the rainbow-ladder approximation uses a bare quark-gluon vertex, 
which violates QCD's Slavnov- Taylor identity (STI). In order to overcome this deficiency, 
physicists are trying their best to go beyond the rainbow-ladder approximation. Much work 
has been done in this direction. Here we just name a few ex amp les: the Ball-Chiu (BC) 
vertex derived from the vector Ward-Takahashi identity (WTI) 18|,ll9|], the CP vertex which 
takes into account some transversal effects 2p| and the vertex derived from the transversal 
WTI 2ll-l23|. etc. As was shown in Ref. j24|, if one deletes the ghost amplitudes and the 
gluon dressing function factor from the STI then the result has the form of a color matrix 
times the WTI structure. Here one notes that the BC vertex ansatz multiplied by the color 
matrix will satisfy such a relation. So in this paper we adopt such a quark-gluon vertex 
ansatz to explore the effect of vertex dressing on the tensor vacuum susceptibility. 
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When one tries to calculate the dressed quark propagator from the DSE using the BC 
vertex, one should construct a consistent kernel approximation corresponding to this ver- 
tex. How to construct systematic and convergent expansions for the kernels of DSE is a 
long-standing unsolved problem. Recently,jm important progress in this problem has been 
achieved in Ref. 



25|. The authors in Ref. [25| have proposed a Bethe-Salpeter kernel which 



is valid for a general quark-gluon vertex. This provides a theoretical foundation for calcu- 
lating the tensor vacuum susceptibility beyond the rainbow-ladder approximation. In the 
present paper we shall use this method to calculate the tensor vacuum susceptibility. 

In order to make this paper self-contained, let us first recall the definition of vacuum 
susceptibilities. In the QCD sum rule external field approach, the QCD vacuum suscep- 
tibility is tightly related to the linear response of the dressed quark propagator coupled 
nonperturbatively to an external current J^{y)Vr{y) = <i{y)^ liyWriy) {qiy) is the quark 



field, r stands for the appropriate combination of Dirac, flavor, co 



or matrices and Vr(?/) is 



the variable external field of interest) Imj. Here following Ref. HI], we adopt the following 
definition for the tensor vacuum susceptibility 



z _ {Tr[a^c^ST ■ ZS] - Tr[ar^(^SoTQ ■ Sq]] 



Z,c(0| : g(OMO) : |0) ' 

where denotes the variable external tensor field, T and Fq ({To)f^i^ = a^y) denote the full 
and free tensor vertex, S and Sq are the full and free quark propagators. (0| : g(0)g(0) : |0) 
denotes the chiral quark condensate. Here it should be noted that Eq. (1) is essentially the 
tensor vacuum polarization, regularized by subtraction of the free vacuum polarization, and 
scaled by the scalar vacuum condensate. It can be obtained by external field differentiation 
of the propagator contracted with a bare vertex. Such a differentiation of a trace of a 
propagator, essentially a condensate, to produce a susceptibility as proportional to the 
associated vacuum polarization has recently been used by some of the same authors in Ref. 



26| for the vector and axial-vector vacuum susceptibility, and also in Refs. 27|, |28|] for the 



scalar and pseudoscalar vacuum susceptibility. 

From Eq. ([T]) we can see that the tensor vacuum susceptibility is closely related to the 
dressed quark propagator and the dressed tensor vertex at zero total momentum. Now we 
turn to the calculation of the dressed quark propagator and the dressed tensor vertex at 
zero total momentum in the DS approach. In the DS approach, the gap equation for the 
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dressed quark propagator S in the chiral limit can be written as 

s{py'^ = Z2ij ■p + 



(2) 



with 

E{p) = Z,j^ g'D,,{p - q)^^^S{q)^ri{q,p), (3) 

where D^y{k) is the dressed gluon propagator and r^(g,p) is the dressed quark-gluon vertex. 
The quark-ghion vertex and quark wave-function renormahzation constants, Zi^2(C^) ^^); 
also depend on the gauge parameter. 

The gap equation's solution has the form 

s{pr' = i^-pA{p\e)+B{p\c) (4) 

and the mass function M{p'^) = B{p'^,('^)/^{p'^X'^) is renormalisation point independent. 
The quark propagator can be obtained from Eq. ([2]) with the following renormalisation 
condition (since QCD is asymptotically free, one can choose this renormahzation condition): 

s{py^\p2=^2 = i7-p- (5) 

The renormalized fully-dressed tensor vertex F^j^ satisfies an inhomogeneous Bethe-Salpeter 
equation: 

.A 

r,,{k, P- C) = Zra,, + / [5(g+)r^,(g, P)5(g_)],,ir[J(g, k- P) . (6) 

Jq 

Here k is the relative and P the total momentum of the quark-antiquark pair; = g±P/2; 
r, s, t, u represent colour and Dirac indices; and K is referred to as the fully-amputated 
quark-antiquark scattering kernel. Zt is the renormalisation constant for the tensor vertex. 

For the specific calculation of i o^is only requires the tensor vertex at P = 0. From 
general Lorentz structure analysis and the asymmetry of the tensor vertex F^j, with respect 
to the indices and z/, we can write down the general form of the tensor vertex 

r^^lP, 0) = a^^^E{p^) + (7^p^ - -f^p^,)F{p^) + Z7 ■ p{'^^Pu - luPp)G{p^) . (7) 

Substituting Eqs. (4) and (7) into Eq. (1), we can obtain the final expression for calcu- 
lating the tensor vacuum susceptibility 

3 /"°° f 
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B{s) 

+P2(s) 



sG{s) 
sA^{s) + B^{s) 



(8) 



where a = (0| : q{0)q{0) : |0) is the two-quark condensate. Here we note that in obtaining 
the above equation, we have made use of the fact that the subtraction term vanishes. 

In phenomenological apphcations, one may proceed by considering the truncation scheme 
for the DSEs and BSEs, especially for the dressed gluon propagator, the dressed quark- 
gluon vertex and the four-point dressed quark-antiquark scattering kernel. The important 
information about the kernel of QCD's gap equation can be phenomenologically drawn by a 
dialogue between DSE studies and results from numerical simulations of lattice-regularized 
QCD 29l432l|. The ansatz that is typically implemented in the quark propagator's gap 
equation can be written as 



Z,g'Dp,{p - q)Tl{q,p) ^ g{{p - qf) D'^^^ip - q)^T^iq,p) , 



(9) 



wherein D^^^'^{i) is the Landau-gauge free gauge-boson propagator, Qi^i"^) is a model effective- 
interaction and rcr{q,p) is a vertex ansatz. 



Over the past 



approximation 



ew years, the most usually used approximation is the rainbow-ladder 



12l4l7l|. where the dressed quark-gluon vertex T^{q,p) is replaced by the bare 



vertex 7^, and in the BS equation the ladder kernel is used. Rainbow-ladder approximation 
is the lowest order truncation scheme for the DSE. It is the nonperturbative symmetry- 
conserving truncation scheme because it satisfies the axial-vector WTI. Models formulated 
using the rainbow-ladder DSE to describe the quark dynamics within hadrons were found 
to provide good and compact descriptions of the light pseudoscalar and vector mesons. 
However, the rainbow-ladder DSE cannot describe well the properties of scalar mesons. So 
physicists are trying to go beyond the rainbow-ladder approximation for years. The key 
points to go beyond the rainbow-ladder approximation are the dressed quark-gluon vertex 
and the four-point quark-antiquark scattering kernel. 

For the dressed quark-gluon vertex, we can employ the BC vertex [l^, [isl 



where 



S^(fc^f) = l[F(fc2) + F(0], 



2 ff2\ 



F{e) - F{e 



(10) 



(11) 



fc2 - £2 ' 

with F = A,B, viz., the scalar functions in Eq. (jlj). Here it should be noted that the BC 
vertex satisfies the vector WTI. 
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Now one should find a kernel consistent with the BC vertex ansatz. This is a difficult 
task that many scientists try to do. Recently, great progress has been done on this aspect. 
The authors in Ref. 25| have found a way to constrain the kernel for the general vertex. 



Following their method, an exact form of the inhomogeneous BSE for the tensor vertex 
^fiu{k, 0) can also be written as 

g^Dap{k - q)—laS{q)T^^{q, Q)S{q)—Tp{q, k) 
+ J^g'D^f^ik - q)^^^S{q)^A^,f,{k,q;0), (12) 

where A^^i3{k,q]0)_js a four-point Schwinger function that is completely defined via the 
quark self-energy [33|, \3^. It satisfies the similar identity as those in Ref. [25 1 



{k - q)i3iA^^p{k, q; 0) = T^^{k, 0) - T^^{q, 0). 



(13) 



Then we can obtain 



iA^^I3{k,q;0) = 2l^[AEiq,k;0) + {-f^l^ - -fJ^)AFiq,k;0)] 

+ {lfiSui3 - juSf,i3)^F{q, k; 0) + 2//37 ■ l{-f^lu - jJ^,)AF{q, k; 0) 
+7 ■ Klf^^ui3 - 'JuSf,i3)T,Giq, k; 0) + 'Jpi'J^.lu - 'yJf,)T.Giq, k; 0) 
+ \Mk^ - g')[7/.(g - k), - - k)^]AG{q, k; 0). 



(14) 



Herein we employ a simplified form of the renormalisation-group-improved effective in- 
teraction proposed in Refs. 
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long-range behavior [s = k"^): 



17l |: viz., we retain only that piece which expresses the 
Q{s) 47r2 



(15) 



This is a finite width representation of the form introduced in Ref. 36|], which has been 
rendered as an 

finite model gap equation. Hence, the regularisation mass-scale can be removed to infinity 
and the renormalisation constants set equal to one. 

The active parameters in Eq. f|T5l) are D and uj but they are not independent. In 
reconsidering a renormalisation-group-improved rainbow-ladder fit to a selection of ground 



state observables [3], Ref. [15|] noted that a change in D is compensated by an alteration 



of UJ. This feature has further been elucidated and exploited in Refs. [l6, 
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FIG. 1: Dressed quark propagator. Left panel - A{p'^), right panel - B{p'^). In both panels, 
Dashed curve : calculated in rainbow-ladder truncation; solid curve : calculated with BC vertex 
ansatz. 

uj G [0.3, 0.5] GeV, with the interaction specified by Eqs. ([9]), ( ITOl) and ( TTSl) . fitted in- vacuum 
low-energy observables are approximately constant along the trajectory 

= (O.SGeV)^ =: m^. (16) 

Herein, we employ uj = 0.5 GeV,D = m^/u = 1.0 GeV^. 

So now with the BC vertex ansatz and the model effective interaction, the equations of 
the DSE for the dressed quark propagator and the BSE for the dressed tensor vertex are 
reduced to a closed system of equations. We can numerically calculate them with iteration 
method. In Fig. [T]we plot the functions obtained through solving the gap equation and in 
Fig. |2] those which describe the dressed tensor vertex. 

It is apparent in Fig. [T]that the vertex Ansatz has a quantitative impact on the magnitude 
and point-wise evolution of the gap equation's solution. That this should be anticipated is 
)lain from Ref. 38|]. Moreover, the pattern of behavior can be understood from Ref. 
39| : the feedback arising through the term in the BC vertex, Eq. ffTOl) . absent in the 
rainbow approximation, always acts to alter the domain upon which v4(p^) and M(p^) differ 
significantly in magnitude from their respective free-particle values. Since -E(p^), F{jp) 
and GijP') are derived quantities, their behavior does not require explanation. We plot the 
integrand in Eq. ([8]) in Fig. |3]for each vertex ansatz. From the figure we can see that there 
is no far-ultraviolet tail in the integrand so that the we do not need regularization here. The 
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FIG. 2: P = scalar vertex, Eq. ([7]): upper left panel - E{p'^), upper right panel - FijP'), lower 
panel - G{p'^). In all panels, Dashed curve : calculated in rainbow-ladder truncation; solid curve : 
calculated with BC vertex ansatz. 

resulting tensor vacuum susceptibilities are 

xlc = 0.05573 GeV-\ xIl = 0.08672 GeV-\ (17) 

The above result shows that the numerical value of the tensor vacuum susceptibility ob- 
tained in the BC vertex approximation is much smaller than that in the rainbow-ladder 
approximation. Here it should be noted that in the above calculations of tensor vacuum 
susceptibility using the effective interaction (15) in the rainbow-ladder truncation and the 
BC vertex, we have chosen the same model parameters for the effective interaction. As is 



shown in Ref. 25|, the amount of chiral symmetry breaking (as measured by the chiral con- 
densate) and related quantities such as the pion decay constant are very different between 
these two truncation schemes. Therefore, when calculating the tensor vacuum susceptibil- 
ity employing the BC vertex, a reasonable approach is to use refitted model parameters in 
the effective interaction (15) in the calculation. Because the active parameters D and u in 
Eq. (15) are not independent, one can refit the model parameters from one physical quan- 
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FIG. 3: Integrand in Eq. ([8]) - Dashed curve : calculated in rainbow-ladder truncation; solid curve : 
calculated with BC vertex ansatz. 



tity, for example, the chiral condensate. Under the BC vertex, the value of the parameter 
D fitted from the chiral condensate is D = | GeV^ (see Ref. 2^]). The results for the 
dressed quark propagator, the scalar functions -E'(p^), FijP'), GijP'), and the integrand in Eq. 
(14), calculated from both the rainbow-ladder truncation and the BC vertex with refitted 
model parameters are shown in Figs. 4 to 6. With refitted parameters in the BC vertex 
approximation, the resulting tensor vacuum susceptibility are 



xlc = 0.07886 GeV 



-1 



xIl = 0.08672 GeV-\ 



(18) 



So, compared with the rainbow-ladder truncation result, the value of the BC vertex 

approximation is reduced by about 10%. Therefore, one can draw the conclusion that in 
the calculation of the tensor vacuum susceptibility in the framework of the DS approach the 
dressing effect of the quark- gluon vertex is not large. 

In Fig. [7] we depict the evolution of the tensor vacuum susceptibility with increasing 
interaction strength, X = Djuj^. The behavior may readily be understood. For X = one 
has a noninteracting theory and the "vacuum" is unperturbed by the external tensor field. 
Hence, the susceptibility is zero. The tensor vacuum susceptibility remains zero until the 
interaction strength X reaches a critical value, X = X^. When X > X^, the tensor vacuum 
susceptibility becomes larger quickly and then goes down slowly for both the rainbow-ladder 
approximation and the BC vertex approximation. Those critical values for the interaction 
strength are: X^^ = 1.93, X^*-^ = 1.41. It can be seen that the critical point in the rainbow- 
ladder approximation is larger than that in the BC vertex approximation. This is easy to 
understand, because the effect of the BC vertex itself amounts to enhancing the interaction 
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FIG. 4: Dressed quark propagator. Left panel - yl(p^), right panel - B{p'^). In both panels, 
Dashed curve : calculated in rainbow-ladder truncation; solid curve : calculated with BC vertex 
ansatz with refitted model parameters. 




FIG. 5: P = scalar vertex, Eq. ([7]): upper left panel - E(p'^), upper right panel - F{p'^), lower 
panel - G{p'^). In all panels. Dashed curve : calculated in rainbow-ladder truncation; solid curve : 
calculated with BC vertex ansatz with refitted model parameters. 
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FIG. 6: Integrand in Eq. 



Dashed curve: calculated in rainbow-ladder truncation; solid curve: 



calculated "with BC vertex ansatz "with refitted model parameters. 




FIG. 7: Dependence of the chiral susceptibility on the interaction strength in Eq. (jl5p : viz., X 
D juj^: dashed curve, RL vertex; solid curve, BC vertex. 



strength. The authors in Ref. [27[ has explained the nature of the critical interaction 
strength which denotes a second-order phase transition. 

For X <Xf., the interaction strength is not sufficient to generate a non-zero scalar term 
in the dressed quark self-energy in the chiral limit. That means below the critical value, 
dynamical chiral symmetry breaking is impossible. The situation changes at X^., for X > X^ 
a B ^ solution is always possible. Moreover, when X < Xc, the interaction strength is also 
not sufficient to generate the non-zero F and G functions in the dressed tensor vertex. That 
is the reason why the tensor vacuum susceptibility remains zero when X < Xc. 

To summarize, using the expression obtained in the QCD sum rule external field approach 
in Ref. we extend the calculation of tensor vacuum susceptibility in the rainbow- 



ladder approximation of the DS approach in Ref. 



11| to that of employing the BC vertex 
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approximation. Here a key problem is how to construct a consistent Bethe-Salpeter kernel 
for a dressed quark-gluon vertex ansatz whose diagrammatic content is unknown. Recently, 
significant progress in this problem was achieved in Ref. 25|. In this paper, following the 



work of Ref. 



251], we construct the kernel for the dressed tensor vertex at P = which 



is needed in the calculation of tensor vacuum susceptibility. Then we perform a consistent 
calculation of the tensor vacuum susceptibility beyond the rainbow-ladder aproximation. 
Our results show that compared with its rainbow-ladder approximation value, the tensor 
vacuum susceptibility in the BC vertex approximation is reduced by about 10%. This shows 
that the dressing effect of the quark-gluon vertex is not large in the calculation of the 
tensor vacuum susceptibility in the framework of the DS approach. In this paper we also 
demonstrate that the tensor vacuum susceptibility can be used to demarcate the domain 
of coupling strength within a theory upon which chiral symmetry is dynamically broken. 
For couplings below the associated critical value and in the absence of confinement, the 
tensor vacuum susceptibility remains zero. This situation changes until the interaction 
strength is larger than a critical point. It is found that the critical point in the rainbow- 
ladder approximation is larger than that in the BC vertex approximation. This is easy to 
understand, because the effect of the BC vertex itself amounts to enhancing the interaction 
strength. 

This work is supported in part by the National Natural Science Foundation of China 
(under Grant Nos. 10775069 and 10935001) and the Research Fund for the Doctoral Program 
of Higher Education (under Grant Nos. 20060284020 and 200802840009). 
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